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INTRODUCTION 

To the Principal Investigator this has been a most satisfying 
contract to have been involved in. The challenge of' studying a 
remote planetary body, of using new data types, of developing 
new analysis procedures has been stimulating and fruitful. At 
the onset of this contract four years ago the Principal Investigator 
felt that his experience studying earth features would be useful 
for interpreting lunar data, ani* that the broadening of one's vista 
would, in turn, be helpful for understanding our home planet, 
the earth. It is gratifying to find that these assumptions have 
proved true. It is also gratifying to find at the conclusion of 
the contract that nearly all the efforts proposed have been success- 
fully accomplished. 

ACKNOWLEDGEMENTS 

I would like to take this opportunity to express my thanks 
to Nat Hardy and Joe Dixon of the JSC TN3/Apollo Photo Data Analysis 
Task Manager Office for encouragement, assistance, and support 
throughout this effort. I am much indebted to W. R. Wollenhaupt of 
the JSC Mission Planning and Analysis Division for completing the 
development of the Extended Sequential Estimation Orbital 
computer program which provides estimates of the acceleration 
components of the space craft from Dopplar tracking data, for 
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providing lunar gravity data processed by this prograun, and for 
helping to explain to a novice the Intricacies of orbital analysis. 
I thank William Sjogren of the Jet Propulsion Laboratory for many 
discussions and for providing preprints and data in advance of 
publication. And, at the Woods Hole Oceanographic Institution, 
Bruce Simon, Leon Gove, Carolyn Dean, Nan Galbraith, and Christine 
Wooding have contributed importantly to the success of this 
contract. Interactions with Wilfred Bryan and Peter Jezek working 
under NAS 9-12564 have also been very ber.eflclal. 

STATUS OF CONTRACT COMMITMENTS 

The status of the principal stages for the investigation are 
summarized below following the underlined stage descriptions from 
the original contract docximent. 

(1) Make gravity attraction calculations at a nominal altitude for 
simple 3-D geometries to simulate craters and maria for initial 
evaluation of effects of topography . 

a . Flat-moon models 

b. Curved-moon models 

A selection of the results from such simple models 
have been published as Figure 1 in Bowin, et al. (1975) . 

These models have proved vesful, for they rather 

quickly showed that mare fill could not match the observed 
gravity measurements unless it is unreasonably thick. 
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(2) Attempt to develop a three-dimensional curved surface 
gravity attraction model program . 

A three-dimensional curved surface gravity attraction 
model progreun has been successfully developed. It is very 
fast computationally and adds negligable computer time to 
that of a flat-surface 3-D program. Our program computes 
both the flat surface and curved surface solutions and 
outputs both values for study. The methodology for this 
progrcwi is given in the Appendix of the paper by Bowin, 
et al. (1975) . 

(3) Evaluate reguirements for processing individual models taking 
into account curvature of lunar surface . 

This evaluation was incorporated into the development 
of the curved surface 3-D gravity computer program, and a 
discussion of some aspects of this matter are discussed in 
the first three paragraphs on page 4949 of Bowin, et al. 
(1975) . 

(4) Make gravity attraction calculations as in (1) , but for 
different altitudes . 

a. Flat-moon models 

b. Curved-moon models 

Evaluate expected vertical gradients to be expected from 
various sized craters and maria. 


Figure 1 of Bowin, et al. (1975) presents gravity 
anomalies computed at both 20 km and 100 km heights above 
the moon's surface, and Figure 4 in the same paper gives 
computed anomalies at all elevations between the surface 
and 300 km height. Consideration of the vertical gradients 
observed and from simple computed models helped lead me 
to the development of the two-body solution for mascon 
anomalies published by Bowin, et al. (1975) . 

(5) Compare vertical gravity gradients produced by homogeneous 
spherical shell lunar mass models with those from model studies 

in (.4).. 

Included in discussion of (4) above. 

(6) Determine anomalous vertical gradients from comparison of 
gravity anomalies obtained at the scime latitude and longitude 
but at different altitudes from the Orbiter and Apollo missions . 

We attempted this, but the rev to rev variations from 
the sequential estimation Apollo radial component data were 
such that this effort was set aside. The low altitude 
sequential estimation Apollo data and the line-of-sight 
Apollo data were received too late in the contract period 
to be incorporated. 
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(7) Determine beat approach for altitude reduction of space- 
craft gravity data . 

The best approach to date Is to subtract a modeled 
acceleration at the position of the spacecraft from the 
observed acceleration of the spacecraft and take the 
difference to be the ” free-air" gravity anomaly at the 
height of the measurement. The large variations in 
vertical gradient over mascon features would introduce 
considerable distortion of the anomalies if a standard 
vertical gradient were applied to all data from many 
different heights. 

(8) Prepare (in conjunction with W. Sjogren and MSC) gravity 
anomaly map of moon's nearside with vertical acceleration values 
reduced to a common datum . 

The difficulties discussed briefly in ( 7 ), the limited 
special distribution of the Apollo data, and the calculation 
of global gravity anomaly maps from orbital rate information 
by Ferrari (1975) lead to the defferement of this effort. 

(9) Cataqorize gravity anomalies with various crater and maria 
paramters . 

Discussed in Bowin (1975) . 

( 10 ) Select specific lunar features for detailed t > *ree-dimensional 
modeling . 

Our analyses concentrated on the Serenitatis mascon 


because it shows evidence of a large mascon, the Apollo 15 
and 17 missions passed nearly over the center of the basin, 
and its simple, approximately circular outline is well 
defined by the surrounding highlands. 

(11) Construct detailed surface and sub-surface three-dimensional 
models utilizing all geological and geophysical data available. 

To be done in conjunction with geologic studies by W, Bryan . 
Compare calculated gravity anomalies with those observed, and 
adjust models in accordance with geologic interpretation and 
gravity differences. This is an iterative, largely empirical . 
process by which limits of likely sub-surface mass variations 
are determined . 

These procedures have led to the development of the 
two-body **olution for mascon anomalies. This solution 
appears to be the geological most reasonable explanation 
for the magnitude of the mascon gravity anomalies. 

(12) Interpret results in terms of crustal structure and origin 
of various crater and maria features on the moon . 

Three papers dealing with such interpretations have 
been prepared under this contract. They are Bowin et al, 
(1975), Bowin (1975), and Bowin (1976). Copies of these 
papers are included in this report. 
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(13) Comparison and interpretation of the similaritlea and 
differences of the gravity fields of the earth, moon and Mars . 

A paper (Bowin, 1976) with essentially this title 
has been submitted for presentation at the VII Lunar 
Science Conference scheduled for March 1976 in Houston, 
Texas. This analysis further suggests an explanation 


for the asteroid belt 



I I- I 


I . J. 
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DATA FILES 

We have prepared a digital data library on magnetic tape 
of single station extended sequential estimation gravity data 
from data received from Johnson Spacecraft Center for Apollo 
missions 14« 15, 16 and 17, These data are tabulated in Tables 
1 and 2, and the format is described in Figure 1. 

PUBLICATIONS 

Three papers have been prepared under support from this 
contract. They are entitled, MASCONS; A TWO-BODY SOLUTION; 
NEGATIVE GRAVITY, ANOMALIES ON THE MOON; and COMPARISON OP 
GRAVITY ANOMALIES TOR EARTH, MARS, AND MOON. The last paper 
is being expanded for submission to a standard journal. 
Reproductions of the three publications follow. 
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IKK - Indicator (-1 for prooont loaar data) 

XSOBC - W.M.O.X. Soareo Coda 

Hat - Say 

RIO - MontA 

Rf* - Yoar 

Km - Hour and Blauto 

SEC - Soeonds 

KAT - XAtitudo in de^raas (aXOO) , 

OLOKO - Lonqituda in de^raoa (aXOO) 

SVCC - Whlela diatance fren center of mm (mS) 

ALTt. - Lasar altitude Ouaa) 

AZ • Atinuth (da^raaal 

SXHC • Inclination (decrees) 

STAC - Tangential acceleration (ingals) 

SHAC - Xoraal acceleration (ngala) 

FA - Xadial acceleration (ngala) 

XKEOH - ntooretical gravity [ngale) 

XSST • Standard deviation of STAC (*100) 

KSSH • Standard deviation of SHAC (*100) 

XSSA - Standard deviation of FA (*100) 

XSSA - Standard deviation of TKSCK (*100) _ 

E1.EV - Elevation of topog r aphy w/ref to radius of 1T3S I0(S 

ELFA - Lasar elevation of topography u/ref to radius of 173S 

ID - Plight and Rev identifier 

ICXS - vehicle identifier 

lATXBY - Sort key of latitude 

LOBCXXY • Sort key of Icngitade 


.deeinal point 


iaplied doeiaal point 


Figure 1. Woo ^3 Hole Oceanographic institution Lunar Gravity Format 
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TABLE 1 


Single Station Extended Sequential 
Estimation Orbital Solutions 


APOLLO 14 


REV 

Number 

of 

Data 

Points 

1 

496 

4 

384 

7 

385 

10 

296 

13 

382 

17 

339 

23 

421 

29 

407 

34 

387 


WHOI 

Source 

Code 

Tape 

Serial 

Number 

File 

Number 

163 

GK14 

1 

163 

6K14 

2 

163 

GK14 

3 

163 

GK14 

4 

163 

GK1.4 

5 

163 

GK14 

6 

163 

GK14 

7 

163 

GK14 

8 

163 

GK14 

9 


II 


TABLE 1. 
(Continued) 


APOLLO 15 

REV 

Nurdaer 

of 

Data 

Points 

WHOI 

Source 

Code 

Tape 

Serial 

Nuihber 

File 

Number 

1 

477 

158 

GK15 

1 

13 

328 

158 

GK15 

2 

16 

422 

158 

GK15 

3 

19 

399 

158 

GK15 

4 

22 

390 

146 

GK15 

5 

25 

422 

146 

GKL5 

6 

27 

414 

146 

GK15 

7 

31 

377 

146 

Gia5 

8 

34 

398 

146 

GK15 

9 

37 

344 

146 

GK15 

10 

40 

309 

146 

GK15 

11 

43 

366 

146 

GK15 

12 

44 

284 

146 

GK15 

13 

54 

362 

146 

GK15 

14 

56 

334 

146 

GKL5 

15 

59 

343 

146 

GK15 

16 

62 

379 

146 

GK15 

17 

65 

299 

146 

GK15 

18 

68 

295 

146 

GK15 

19 

70 

348 

146 

GK15 

20 

71 

318 

146 

GK15 

21 

72 

371 

146 

GK15 

22 
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TABLE 1. 
(Continued) 
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J _L 


APOLLO 16 

Number 


REV 

of 

Data 

Points 

WHOI 

Source 

Code 

Tape 

Serial 

Number 

File 

Number 

1 

480 

160 

6K16 

1 

3 

354 

163 

6K16 

2 

6 

368 

163 

GK16 

3 

9 

365 

163 

GK16 

4 

12 

346 

163 

GK16 

5 

17 

429 

160 

GK16 

6 

22 

329 

160 

GK16 

7 

27 

425 

160 

GK16 

8 

32 

395 

160 

GK16 

9 

37 

213 

160 

GK16 

10 

42 

399 

160 

GK16 

11 

47 

431 

160 

GK16 

12 

52 

430 

160 

GK16 

13 

57 

400 

160 

GK16 

14 

61 

422 

160 

GK16 

15 


APOLLO 17 

TABLE 1. 
(Continued) 

Number 

of WHOI 

Data Source 

Tape 

Serial 

File 

REV 

Points 

Code 

Number 

Number 

1 

412 

158 

6K17 

1 

3 

379 

165 

6K17 

2 

5 

384 

165 

GK17 

3 

6 

389 

165 

GK17 

4 

7 

383 

165 

GK17 

5 

8 

375 

165 

GK17 

6 

9 

384 

165 

GK17 

7 

13 

396 

159 

GK17 

8 

14 

220 

143 

GK17 

Q 

15 

434 

143 

GK17 

10 

21 

200 

160 

GK17 

11 

23 

216 

158 

CI'.17 

12 

24 

249 

143 

GK17 

13 

27 

448 

14J 

GK17 

14 

28 

55 > 

143 

GKl 7 

15 

33 

200 

158 

GK17 

16 

38 

539 

143 

GK17 

17 

42 

392 

158 

GK17 

18 

47 

210 

159 

GK17 

19 

49 

310 

143 

GK17 

20 

51 

395 

160 

GK17 

21 

59 

415 

158 

GK17 

22 

62 

523 

143 

GK17 

23 

64 

441 

158 

GK17 

24 

65 

380 

143 

GK17 

25 

66 

475 

143 

GK17 

26 

71 

602 

143 

GK17 

27 

73 

443 

158 

GK17 

28 

74 

492 

143 

GK17 

29 


1 

1 

^ .. 

TABLE 

1 . - 4 . — ■ 

2. 


1 



Single Station Solutions With 
Laser Altitude 



APOLLO 

MISSION 

REV 

Number 

of 

Data 

Points 

Number of 
Data Points 
with Laser 
Altitude 

WHOI 

Serial 

Number 

File 

Number 

15 

16 

422 

121 

GK15 

3 

15 

22 

390 

227 

GK15 

5 

16 

17 

429 

51 

GK16 

6 

16 

37 

213 

54 

GK16 

10 

16 

47 

431 

118 

GK16 

12 

17 

14 

220 

61 

GK17 

9 

17 

15 

434 

118 

GK17 

10 

17 

24 

249 

75 

GK17 

13 

17 

27 

448 

79 

GK17 

14 

17 

28 

555 

161 

GK17 

15 

17 

38 

539 

160 

GK17 

17 

17 

49 

310 

107 

GK17 

20 

17 

62 

523 

147 

GK17 

23 

17 

65 

380 

44 

GK17 

25 

17 

66 

475 

106 

GK17 

26 

x7 

71 

602 

193 

GK17 

27 

17 

74 

492 

119 

GK17 

29 




15 


i 


1 





VOL. 80, NO. 35 




A 


A 



JOURNAL OF GEOPHYSICAL RESEARCH 


Mascons: A Two-Body Solution 

Carl Bowin and Bruce Simon 

Woods Hole Oceanographic Inslilulion, Woods Hole, Massachusetts 02343 

VV. R. WOLLENHAUPT 
NASA Johnson Space Center, Houston. Texas 77083 
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Almost all of the muss distributions ihut have been proposed to account for Ihe large positive gravity 
anomalies associated with lunar mascons have assumed single-body sources of a mass excess. In the case 
of mare fill with a reasonable density contrast (-1-0.5 g/cm*) with crustal material, a fill thickness of about 
16 km for Mare Serenitatis i< thus required to account for the observed gravity values at lOO-km height. 
Such a great thickness would require a 16-km-deep hole prior to filling, and such a topographic depression 
is inconsistent with the depths of the topography of the Mare Necturis and Mare Oriental basins, which 
have little fill, and with estimates of mare thicknesses based on buried crater dimensions. A two-body 
mascun solution, however, requires only about a 2-km thickness of fill and a 12-km rise of a lunar Moho 
beneath Mure Serenitatis to account for observed gravity anomalies. The mantle dome results from 
an uprising of mantle material beneath the mare basin, bringing the impact crater to near isostatic 
equilibrium Two kilometers of fill is inferred to have accumulated later, when the crust became rigid 
enough to sustain the loud. Together the fill and the dome account (at about 20% and 80%, respectively) 
fur the magnitude of the observed mascon anomalies. This type of two-body solution cun account for 
greater magnitude mascon gravity anomalies by proportional increases in the fill thickness. The top 
of the mantle dome or plug is placed at 60-km depth to match observed seismic velocity structure. This 
mascon structure has an anomalous gravity field that is in agreement with the maximum magnitude of 
anomalies observed at several heights above Mare Serenitatis. Hie thickness of fill would be greater if 
the basin floor had subsided under the load of early fill material. 


Introduction 

The origin of mascons (lunar mass concentrations) has been 
of interest since their discovery by Muller and Sjogren (1968). 
Almost all of the mass distributions that have been proposed 
to account for the large positive gravity anomalies associated 
with mascons have been isostalically uncompensated struc- 
tures. Iixa:'i,'les of such structures are buried iron meteorites 
[Stipe, I968J, >he pooling of a melt from an impacted meteorite 
[Urey, 1968], a mantle plug [Wise and Yales, 1970], i • the 
mare fill ilself [Conel and Holsirotn, Baldwin. 1968; I'hil- 
lips et al., 1972: Wood. 1970; Booker el ai, 1970], 

These models have all been single-body sources of a muss 
excess to account for the mascon gravity anomaly. Phillips el 
al. (1972, p. 7113] briefly discussed the po.ssibility of a two- 
body mass anomaly consisting of both marc fill and a mantle 
plug. However, they did no! consider that a plug would con- 
tribute significantly to the gravity etfect observed. Arkani- 
Named [1973] propo.sed development of a 20-km-high mantle 
plug beneath Mare Imbrium immediately following impact, 
and he examined temperature variations with depth for such a 
structure. The gravitational attraction of this structure was not 
calculated. 

Only Nulme [1972] has computed the gravitational attrac- 
tion of a two-body mass anomaly distribution, and it is unique 
in being the only mascon struc. are so far proposed that could 
be in isosiatic equilibrium (Figure \f). The specific structure 
that he presents, however, is not in perfect isostatic equilib- 
rium (\V. L. Sjogren, personal communication. 1974). because 
the crustal layer is thinner beneath the center of the mare and 
thicker at the margin than is required for equilibrium. 

Wise and Yates [1970] proposed that a mantle plug formed 
following the formation of u deep crater by meteorite impact. 


The rise of this plug and the concomitant volcanic fill reestab- 
lished isostatic equilibrium but still left a crater some 6.4 km 
deep. Because a quantity of volcanic material is less dense 
when it is molten than when it is solid. Wise and Yates postu- 
late that continued flooding of the crater by lava occurred until 
the excess of lithostatic pressure over hydrostatic pressure 
could lift the magma no further. The continued flooding pro- 
duced a supe. iv 0 .static fill. In their analysis Wise and Yales 
inferred that the volcanic fill has the same density as the nor- 
mal lunar crust, and hence they concluded that after the mare 
basin filled with this material, . . the fill would disappear 
gravimetricully. and the residual positive anomalies in the 
gravity picture would be due to the Moho plug plus any topo- 
graphic effects.' This is the case because in their model the 
density of the .superisostatic fill is the same as the density of the 
.surrounding cru.st. By assuming a density contrast of 0.5 g/cm* 
between the plug and the lunar crust the gravitational attrac- 
tion of several configurations for the mantle plug was com- 
puted to match the Lunar Orbiter data for Mare Imbrium 
[Muller and Sjogren, 1969]. For the model assuming a crustal 
thickness of 61 km the lieight of the mantle plug is 31 km, and 
the depth to the lop of the plug is 30 km. 

In this paper we present theoretical free-air gravity anoma- 
lies that occur over a variety of simple structural con- 
figurations for reference and us an aid to the understanding of 
the constraints imposed by the observational gravity data. 
These calculations take into account the curvature of the lunar 
surface. The computer program which performs these com- 
putations is described in the appendix. By means of the simple 
structures presented, available observational gravity data, and 
information on the depth of mare fill and rock densities we 
argue that a single-body mass anomaly is not adequate to 
explain the mascon gravity anomalies, and we present a pre- 
ferred two-body solution. 


Copyright ® 1975 by the American (ieophysical Union 
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Hg. I. CruM sections of mare basin structures and their gravity anomalies. Structures on the right are in isostatic 
equilibrium, and those on the left are not. Lunar crust is indicated by stippled pattern (assumed density of 2.9 g/cm’). Mare 
fill and source of mascon anomaly are indicated in black (assumed density of 3.4g/em’). The unfilled mare basin (density 
of zero) is left unpatterned, as is the lunar mantle (density 3.4 g/cm’) in part/. In parts g. c, and rf the thickness of the fill is 
8.53 km. and the overlying depressions are 1.47, 0.67. and 2.97 km deep, respectively. 


SiMHi.E Mare Models 

The frec-uir unumnlies that would be associated with several 
simple mare basin configurations are shown in l-igure I. The 
configurations in the right-hand column are isostatically com- 
pensated. and those in the left-hand column arc not. The free- 
air anomaly profiles are not diagrammatic but are accurately 
calculated from the models by taking into account the curva- 
ture of the moon by means of the computer program described 
in the appendix, hree-air anomaly profiles (Figure I ) aie given 
at both 20- and lOO-km height in order to give the reader a 
better feeling for the change of the anomaly magnitude with 
height. 

A hypothetical sequence of marc evolution follows to illus- 
trate changes in the gravity field that would occur over a 10- 


km-deep basin. This hypothetical sequence does not involve 
any deep anomalous mass distributions; it involves only those 
of the basin and its fill. Prior to the filling of the basin (Figure 
Ifl). large negative gravity anomalies would occur. As the 
basin began to fill, the magnitude of the negative anomaly 
would decrease, and a local gravity high would develop in the 
center (as occurs in Figure \d but with much less than 10 km 
of fill). At a certain stage of filling (F'gure Ig) the thickness of 
the fill would be such that its mass in excess of that of the 
equivalent volume of crust would just equal the mass deficiency 
of the remaining void above the fill, and the structure would 
be in iso.stulic equilibrium. Additional filling would lead to an 
excess of mass and to positive gravity anomalies (Figure If). 
Upon complete filling of the basin (Figure Ifc) the grea’.. • 
positive anomaly would occur. We might now imagine that 
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the basin floor would begin to link. This would lead to a 
reduction in the magnitude of the positive gravity anomaly 
(Figure Ir). Sinking might continue until isostatic equilibrium 
again occurred (Figure Ig). If sinking continued further, nega- 
tive gravity anomalies would become more and more pro- 
nounced (Figure \d). Of course, all these stages need not be 
imagined, and some intermediate stages might be bypassed 
in other hypothetical sequences. The models of Figure I are 
illustrative of the free-air anomalies that would be expected 
for a variety of possible evolutionary paths. 

For the calculation of gravity anomalies for the models of 
Figure I there are two general approaches. One can use the 
actual assumed densities of each of the portions of the model, 
sum their attractions, obtain the difference between the total 
attraction and the known or assumed anomaly value at some 
location, and add that difference algebraically to the total 
attraction value at other locations to obtain anomalies. Or one 
can assign only the value of density contrasts to each portion 
of the model. The two methods provide identical results except 
in the case in which the strong gravity gradients that occur at 
the truncated edges of the model also have an appreciable 
effect on the portion of interest in the model. This condition is 
the case for three-dimensional models of circular maria on the 
moon’s surface. Accordingly, better results are obtained if 
density contrasts are used in the computations. The density 
confasts are those mass distributions that differ from tbe 
mass tribution of a moon composed of homogeneous con- 
centric $, ‘'erical shells. 

For all the gravity calculations for all the models of Figure I 
the gravitational attraction of the mare fill is only that of the 
density contrast between the fill and the crust material. For the 
whole rock density to be used as a density contrast in com- 
puting its gravitational attraction the fill would have to be 
isolated laterally from the surrounding crust. That is, it would 
have to be above the normal level of the lunar crust. This type 
of model is not shown in Figure I because the mascon anoma- 
lies occur over depressions, rather than highs, in the lunar 
topography. 

In this paper we assume that the normal level of the lunar 
crust is equal to that of the mare surface. This, of course, may 
not be true. The mure surfaces may have sunk, or the mare 
depression may not have been filled as high as the normal level 
of the lunar crust. If the latter situation is. assumed, then a 
volume of space having a density of zero (or a contrast of 
minus the lunar crust density) should be incorporated into a 
model used for gravity computations. Zero density leads to 
very large mass deficiencies, as is readily clear from the large 
negative free-air anomaly of Figure lo. Thus to produce the 
mascon anomaly, the subterranean mass excess must be 
greater than the value obtained if the normal level is taken at 
the mare surface. Because the observed gravity anomalies (Fig- 
ure la) do not indicate the occurrence of large positive values 
over the highlands, we assume that the muss of the topo- 
graphic elevation of the highlands above the mure surface 
generally is largely compensated at depth by a mass deficiency. 
This is probably due lo a thicker Ci ast, which for the purposes 
of this paper has not been included in the structure models. 

Lunar Gravity Data 

Mascons w ere discovered by Muller and Sjogren (196k) dur- 
ing their processing of Lunar Orbiter tracking data. Several 
recent studies of lunar mascons have concentrated on an analy- 
sis of the gravity field over Marc Screnitatis [Conel and Hoi- 
slroni, 1968; lluhne, 1972; Sjogren el al., 1972; Phillips el ai. 


1972]. Serenilatis is of interest for several reasons. It shows 
evidence of a large mascon, the Apollo 15 and 17 missions 
passed nearly over the center of the basin, and its simple, 
approximately circular outline is well defined by the surround- 
ing highlands. 

Mullerelal. (1974) determinedu maximum free-air anomaly 
value of -1-230 mCal at a height of 15 km over Mare Sereni- 
tatis. This anomaly was based on performing spline fits to the 
converged least squares Doppler residuals in determining line- 
of-sight acceleration [Sjogren el ai. 1974). At the longitude of 
Serenitatis the line-of-sight acceleration should be only about 
5.5% lower than the moon's radial component of gravity accel- 
eration. 

The g'u>':\y data for this analysis were obtained with a 
relatively new computer program developed at the Johnson 
Space Center. This program uses extended sequential estima- 
tion technique for processing the Doppler data to solve for 
the spacecraft state vector. During computation the spacecraft 
state vector and the components of the unmodelcd accelera- 
tions are simultaneously determined. The spacecraft stale vec- 
tor consists of the components of the position and velocity 
vector of the spacecrai' at each data measurement lime. Dop- 
pler residuals are obtaii'.ed by subtracting from the observed 
Doppler signal a mcJclcd acceleration based on a gravity 
mode) known as the LI model. This model contains the follow- 
ing spherical harmonic coefficients: J„ - 2.07108 X 10 *, C« 
= 0.20716 X 10-*,^*.= -0.21 X I0-‘,C„ -=0.34X 10 *. and 
Cu * 0.02583 X 10'*. The unmudeled accelerations arc deter- 
mined for three components: radial (toward the center of the 
moon), tangential (in the direction of the orbit path and nor- 
mal to the radial component), and normal (perpendicular to 
the orbit path and at right angles to the other two com- 
ponents). The radial acceleration is equivalent to free-air 
anomalies at the elevation of the observations. 

The resolution of the program is primarily limited by the 
noise on the trucking data. The major source of this noise 
appears to be short-term differences between the cesium fre- 
quency standards at the individual tracking stations. When 
one is simultaneously processing multiple-station data, the 
noise in the estimated accelerations is of the order of 15-20 
mGal (Figure 26). If data from only one tracking station are 
processed, the noise appears to be of the order of 5 mGal 
(Figure 2a). The revolution to revolution consistency is nor- 
mally within about 5 mGal. Occasionally, the difference is up 
to 20 mGal. The loss of parallax information in the single- 
station solutions results in attenuation of the estimated or 
recovered accelerations. This is particularly true for the limb 
regions, say, >60® longitude. For Screnitatis the attenuation 
was only about 15 mGal. Since this article is primarily con- 
cerned with the central region, it was decided to use the single- 
station solutions. It is estimated that the recovered accelera- 
tions in this region, derived from single-station solutions, arc 
good to 40-60 mGal for the individual components and to 
about 20-30 mGal for the total recovered acceleration. The 
unmodeled radial accelerations estimated by this program are 
consistent with the Jet Propulsion Laboratory (JPL) line-of- 
sight accelerations within the limitations of each technical ap- 
proach. 

l.arge negative free-air anomalies are observed in the Apollo 
Command Ship Module (CS.M) tracking data adjacent to the 
positive anomalies associated with mascons. The negative 
anomaly is especially pronounced between MarcCrisium and 
Mare Serenitatis. The origin of these lows is not known. They 
appear to be real features and not artifacts of the processing 
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Fig. 2. I- m-uir gravity anomalies from Apollo 13 revoli lion 31. (a) Single-tracking station solution, (h) Multiple-tracking 

stat on solution. 


techniques, because u large negative free-air anomaly ( -21.3 ± K 
mGal) was derived [Sjogren el ai. 1974] at the Apollo 17 
landing site from the Lunar Traverse Gravimeter results \Tal’ 
want el ai, 1973) and it coincides in location with the negative 
anomaly minimum found in the orbit tracking data (Figure 3). 

The negative gravity anomalies adjacent to mascons have 
commonly been inferred to ring the mascon completely [IVise 
and Yaie.\, 1970; Uulnie, 1972; Arkani-Uamed, 1973). This as- 
sumption docs not appear to be generally true. F igure 3 shows 
the location of free-air anomaly minima on the sides of Marc 
Serenitatis as observed in the Apollo 15 revolution having a 
height of more than !U0 km. This figure is based on single- 
station solutions, for which the crest and trough of anomalies 
arc more sharply defined than those for multiple-station solu- 
tions. These data clearly show that the free-air minimum on 
the east side of Serenitatis Joes not continue around to the 
southern margin of the marc. Alsu, the positive anomaly of the 
Serenitatis mascon extends beyond the southern margin of the 


marc (although it is much diminished in magnitude there) at 
least as far as the distance from the though of the minimum to 
the cast mtirgin. A map of line-of-sight accelerations over 
Marc Serenitatis and the region to the south and east prepared 
by W. L. Sjogren (personal communication. 1975) from 
Apollo 15 eSM data for revolutions 15-70 is consistent with 
the above interpretation. 

Mark Fit i. Thu knkss am> 

Ol NSttV CoMRASt 

Gravity anomalies result from a distribution of mass anoma- 
lies and hence depend upon the volume, density, and location 
of the mass anomalies. Thus from gravity data alone it is 
impossible to define uniquely the thickness of the marc fill and 
its density contrast with the crustal material. That is why the 
inverse solutions presented by Phillipf el ai (1972. Figures 2 
and 4) are given in terms of surface densities for various 
depths. For zero depth a surface density of about K (g/cm*) km 
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I ig. 3. I.uculiuii uf rr«i;-air unomal> minima and maviimim in Marc Screniiati> region. Minima arc indicated h> d.i%hcs, 
and maximum b) pluscx. D.itu houndaric% from lliotc Apullo 15 orliii revolutions 1-73 having more th.in lOU-km 
height arc shown h) thin solid lines. 


u us determined for the Screnitulis mascon. I'liis value curre- Cunel timl HoLurom ll%8| matched the pcak-to-trough 
spiuids, for example, to a till thickness of It km for a densit) amplitude (from about l I7U to - KO mCJal) for the Screnitatis 

contrast of 1.0 g/cni’ or to a thickness of 16 km for a densit) mascon and determined a (ill thickness of 14 km for a densit) 

contrast of 0.5 p/cm*. contrast of l.l g/cm’. They recopni/e that this is an unusually 
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high density coiitrusi lo accept. They also determined u till 
thickne!» of 30.8 km for a density contrast of 0.5 g/cm*. Using 
an infinite slab approximation, Wood [1970] calculated that a 
I'km-thick layer of density 3.3 g/cm’ would give rise lo a 
gravity anomaly of about +140 mCal and hence would be 
comparable to mascon anomalies. As was discussed pre- 
viously, the density contrast, not the rock density, must be 
used in such calculations. In addition. Wood's estimate of the 
gravitational anomaly is too high because it was bused on 
the assumption that the l-km-thick layer is an infinite slab. 
The gravitational attraction of an infinite slab is independent 
of the distance from the surface of the slab. Taking into 
account the finite dimensions of a l-km-thick layer (of 
3.3-g/cm’ density) of mure fill of 320-km radius, we calculate 
that a gravitational anomaly of 98 mCal at lOU-km height 
would be produced, a reduction of 30% from the infinite 
slab approximation. 

From the preceding review it is clear that the fill in a mure 
basin such as Screnitatis would have to be about 16 km thick if 
it hud a density contrast of 0.5 g/cm* with the normal crustal 
material in order to account for the mascon gravity anomaly. 
Such a great prefill basin topography seems highly unlikely. 
This thickness would be reduced to about 8 km if the density 
contrast were 1 .0 g/cm’; however, such a high density contrast 
is unlikely. 

Density measurements on lunar samples are unfortunately 
very limited. A summary of published and unpublished values 
is given by Tahvani el al. |l973j. If one excludes those measure- 
ments for which the volume was determined from aluminum 
foil models, the densities of the basalt samples range from 2.8 
to 3.4 g/cm’, and those of the breccia samples range from 2.2 
to 3.0 g/cm*. The three basalt samples measured by hydro- 
static weights all provide a dcn.sity of 3.4 g/cm’. Those samples 
for which porosity was also determined indicate that intrinsic 
densities range from 3.25 to 3.49 g/cm’ for the three basalts 
and from 2.99 to 3.14 g/cm’ for the live breccias. Although 
near-surface rocks may have low bulk densities because of 
iheir vugs, vesicles, or fractures, within a few kilometers of 
depth they should nearly approach their intrinsic grain den- 
sity. Thus for constructing structure models for gravity studies 
the assumption of a basalt density of 3.4 g/cm’ and a density 
contrast of 0.5 g/cm’ with highland breccia material (rock 
density of 2.9 g/cm’) is supported by the available measure- 
ments. The density contrast of 0.5 g/cm’ may be slightly high. 

For a density contrast of 0.5 g/cm’, if the thickness of mare 
fill in Serenitatis is less than 16 km, the icsulting maximum 
free-air anomaly over the mare will be less than that observed 
(Figure \h). Since the maximum anomaly value for individual 
mascons appears lo be proportional to the si/e of the marc 
basin, the last statement is probably true for the other mascons 
as well us for the one at Serenitatis. 

Evidence from crater counts and from radioactive age dat- 
ing of returned samples suggests that considerable time (per- 
haps about 200 X 10* years) elapsed between the formation of 
the large maie basins and their filling by mare basalts [Bald- 
win, 1963. \91\ . Uimekeet al., 1973). The existence of mascons 
appears lo provide evidence that the rigidity of the lunar cru.s| 
was increasing during the same interval [Baldwin. 1971; Phil- 
lips el al.. 1974; Muehiberger, 1974). A lO-km-dcep com- 
pensated hole in the lunar surface would have a free-air anom- 
aly of about -3S5 mfial (Figure \e). and a 16-km hole would 
have an even greater negative gravity anomaly. Since gravity 
anom.dies that large arc not observed on the moon and since 
the present lithosphere is considered to be more rigid now than 


when the marc basins were formed, the premare fill depression 
was probably only a few kilometers deep at most. This deduc- 
tion is compatible with estimates of the thickness of mure fill in 
shelf regions on the basis of inferred dimensions of buried 
craters [Marshall. 1963; Baldwin, \91\', Oe Hon, 1974). These 
results suggest thicknesses of mure fill that are less than 2 km. 
Probably less than I km of fill iKcurs in the shelf area of Mare 
Imbrium [Baldwin. 1971). The average thickness within Mare 
Tranquillitatis is 5(X)-600 m, with a maximum accumulation in 
excess of 1200 m in the vicinity of Lamont crater [De Hon, 
1974). The flooded inner basin of Necturis has a fill that is at 
least 1.2 km thick. The maximum thickness is undetermined 
but may not exceed 1.5- 1.6 km [De Hon, 1974). Mare Nec- 
taris and Mare Oriental have little fill in comparison with 
Crisium, Serenitatis, and Imbrium but do not show a pro- 
nounced increase in depth as the center of the basin is 
approached. From this comparison and from the estimates 
of fill thickness from buried craters it is very unlikely that 
great thicknesses of marc fill occur in the centers of the 
mare basins. 

The small mascon gravity high over Tranquillitatis corre- 
lates in location with the thicker section of fill material at 
Lamont [Phillips and Saunders, 1974). Again, however, the 1.2- 
km-thick section at Lamont can only yield the mascon gravity 
anomaly value if a high density contrast (of about 3.15 g/cm’) 
is assigned to it [Phillips and Saunders, 1974). Such a density 
contrast is not realistic, and even though the nreftll topography 
may have been in isostatic equilibrium, it would not be correct, 
as was explained previously, to use the rock density of the 
superisostatic fill in gravity model calculations. Mence the fill 
alone cannot explain the mascon gravity anomaly. 

From the foregoing discussion it seems to be amply demon- 
strated that the mascon gravity anomalies are not produced 
simply by the occurrence of a superisostatic fill in mare basins. 
The fill alone is too thin to account for the large positive 
anomalies observed. It is also indicated that the mass distribu- 
tion causing the mascon anomalies must be out of isostatic 
equilibrium. 

Two-Booy Mascon Soi.ltion 

A simple solution to account for the mascon anomaly can be 
achieved by considering the mass excess to occur in two bodies 
rather than in a single body. This solution conceptually is 
much like that of Wise and Yales |1970], except that here the 
fill, as well as the mantle plug or dome, is considered to con- 
tribute to the mascon anomaly. Consider the premare basin to 
have had a structure like that in Figure Ir, only not nearly so 
deep. This structure would be in isoslatic equilibrium because 
a mantle dome having excess mass rose to balance the mass 
deficiency of the depression in the crust surface caused by a 
large impact. For the model of Figure 4, 14 km of crust is 
inferred to have been blasted away by the impact, an action 
resulting in a thinned crust. Note that although this structure 
is in isoslatic equilibrium, the free-air anomaly would not he 
zero over the crater basin. Instead, negative values would 
occur at the center. The magnitude of the minimum depends 
considerably upon the thickness of the crust. As the mantle 
dome approaches the .surface more closely , the magnitude of 
the minimum diminishes; if the dome lies at greater depths, 
the minimum increases. This is known as the topographic 
effect, which rc.ults because gravitatiorral attraction varies 
inversely as the square of distance, and the deficiency of the 
crater is closer to the gravity observation site than is the 
excess of the mantle dome. 
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Fif. 4. Two-body muKon itructurc model. Figure fhowrt a crou icction through the center of a radially lymmetrical 
three-dimensional structure. Isounomaly lines show computed Tree-air anomalies above this structure. Numbers in the 
model are densities in grams per cubic ceniimetcr. Shaded bodies are the mass excess producing the maKon gravity 
anomaly. Values in space are observed Tree-air anomalies Trom Apollo 15 revolutions I (higher) and 22 (lower). 


It is likely that the rigidity of the lunar crust increased with 
time, following the formation of the mure basin. Thus at the 
time of the surface eruption of the mare fill volcanism the 
rigidity of the crust would be such that it could sustain the loud 
due to the fill without the occurrence of major isostatic adjust- 
ments at depth. A structure model for Mare Serenitatis based 
on this evolutionary outline is shown in Figure 4. Prior to 
filling, the basin was 2 km deep in crustal material having a 
density of 2.9 g/cm* and is inferred to have been isostatically 
compensated by a 12-km rise of lunar mantle material having a 
density of 3.4 g/cm*. This depression would have had a min- 
imum free-air anomaly value of about -44 mGal over the 
center at lOO-km height. During the interval 3.9-2. 1 X lU* 
years ago. filling of this basin by flooding of basalt lava oc- 
curred [Huneke el al., 1973]. Although there is evidence for 
some subsidence of Mare Serenitatis [Bryan, 1973; Muehiber- 
ger, 1974). the crust is assumed to have been sufficiently rigid 
to support the main load represented by this volcanic fill. The 
magnitude of the excess loud is CKO kg/cm* beneath the center 
of the basin, when a lava density of 3.4 g/cm’ is assumed. 
Mowever. its anomalous gravitational attraction results only 
from the density contrast of the fill with the surrounding crust 
(taken to be -f-0.5 g/cm* in the model of Figure 4). Since the 
large mass deficiency that previously existed because of the 
unfilled depression has now disappeared, the gravitational at- 
traction of the density contrast between the mantle plug and 
the adjacent crust is no longer canceled by the negative attrac- 
tion of the topographic depression. Hence the plug's anoma- 


lous attraction adds to that of the fill, the combination result- 
ing in a large positive gravity anomaly over the mare. 

At a height of 17 km above the center of the model (Figure 
4) the fill contributes 21.3% of the total anomalous gravita- 
tional attraction, and the mantle dome 78.7%. The percentage 
contribution of the fill diminishes gradually with increasing 
height, so that at 300 km it is 18.9% and that of the mantle 
dome is 81.1%. If a deeper (or shallower) original mare basin 
were assumed, then a greater (or lesser) total anomalous gravi- 
tational attraction would result from the model, but the per- 
centage contributions of the mantle dome and fill would nut 
change greatly. Differences in the assumed density contrasts 
could be accommodated within limits by different original 
mure basin depths and yet could have nearly the same total 
anomalous gravitational attraction for the models. The fill and 
the mantle dome could also be assumed to have different 
densities. As greater confidence in the precision accuracy 
of the lunar gravity anomalies is obtained, it nijy become 
possible to limit more the relative contributions of shallow and 
deep sources for the mascun anomalies. Without the gravita- 
tional contribution of the fill the mantle dome alone would nut 
match the observed values. To match the total anomaly of 
Figure 4 would require that the mantle have a density contrast 
with the crust of 0.6.) g/cm*. But of course if it did. then it 
would be smaller in volume and hence would have less gravita- 
tional attraction, since the dome would not need to rise as high 
in order to provide compensation for the mass deficiency of 
the basin. It therefore seems most likely that both mantle 
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dome und (ill conirihulions ure required In otpbin muwon 
fruviiy anomulic*. 

In the model of I'igurc 4 (he lop of the munlle dome mus 
locuied al 6U*km depth to be compalihle with the limited 
information aviiilaMc on the >ei»mic velocity structure nl the 
nuHin [iaminlein vt al., I974|. No marked change in Kivmic 
cttmpreMional velocity is noted at the I* to J-km depth corre* 
sptmding to the base of the mare fill. We have no esplanutinn 
for this possible disagreement with our inascon model (Figure 
4). except to suggest that perhaps the fractured nature nf the 
shallow material may mask the velocity contrast that would 
be expected nl the assumed density contrast. 

The two-body mascon model shown (Figure 4) is in agree- 
ment with the large positive free-air gravity anomaly values 
determined from lunar-orbiting spacecraft. Single-station 
solution values for Apollo 15 revolutions I and 22 are given 
in Figure 4. I'he theoretical anomaly values are on a plane 
across the center of the circular symmetrical structure model. 
The observational data are plotted with respect to the maxi- 
mum anomaly value along its particular orbit. Revolution I 
passed closer to the presumed center of the Mare Sereniiatis 
masctm than did revolution 22. However, the data for these 
revolutions arc considered to be representative and to indicate 
the general agreement between the two-bsxiy maKon solution 
und the observational data. I*reliminary processing of low- 
altitude (about 2U km) revolutions for Apollo 15 (revolutions 
3-12) indicates agreement with the theoretical values of the 
model (Figure 4). 1 he agreement between the theoretical und 
the observational data diminishes toward the periphery of the 
marc basin because the model docs not account for the large 
negative anomalies that Dank the basin along the trajectories 
of the Apollo 1 5 und 1 7 missions. 

The loading of this mascon structure model (Figure 4) is 6S0 
kg/cm'. 1 his value is slightly less than the 800-kg/cm' loading 
value derived by Sjogren cial. (1972) for Mure Nectaris and by 
Phillips el al. (1972) for Mure Serenitutis. This lower value 
arises because the radial component derived from the sequen- 
tial estimation program generally yields lower-magnitude 
anomalies than does the JPL spline fit program fur the total 
acceleration. Increasing the fill thickness and the elevation of 
the mantle dome of our model (Figure 4) can make it com- 
patible with higher mascon anomaly values. The percentage 
increase in fill thickness is approximately the same us the 
percentage increase in the anomaly. For example, a 150-mCial 
m.iscon anomaly al KlO-km height would be produced by a fill 
thickness of 2.9 km. 

In addition to possible errors in assuming ruck densities or 
matching calculated anomalies to a gravity value tha' is too 
low, the lilt may actually be greater than 2 km thick for an- 
other reason. If sinking of the basin flavor occurred in the early 
stages of basin tilling in response to loading due to the fill, then 
some portion or all of the curly lill. dependtng upon the 
amount of sinking, would be compensated by the depression 
of the lltMir. During at least the later stages of marc filling the 
lunar crust became rigid enough to support the mass of the 
new fill without sinking. If all the filling occurred after this 
condition had been .u hieved, then only about 2 km of lill is 
required (together with the mantle dome) to account lor the 
mascon gravity anomaly. 

I he structure model of I igurc 4 provides a two-body mas- 
con solution ih.it can account for the gravity field - ' ''ar ob- 
served at several heights over the central part of M ' eni- 
latis I he tliickiicss of the mare lill and the density values 
assumed are consistent with the data available. 


AmNinx: CiRAViTATinNAi. Attiia( tion Fr<»vi a 
T iiRi-.H-DisihNsioNAi. B(N)Y on a Cuiivh) Si;iii'Ace 

As the si/e of a body whose gravitational anomaly one 
wishes to compute becomes relatively large in comparison with 
the radius of curvature of the sphere in which it resides or if the 
anomaly is desired at points distant along the surface, it be- 
comes necessary to employ equations which account for this 
curvature. Following the procedure of TaUanI and Ewing 
1 1960] the body is divided into polygonal contours, which are 
subdivided into triangular segments. As the attraction from 
each triangular lamina segment is exactly computed, a curved 
correction term is added to it, producing a net anomaly for a 
curved triangular lamina segment, '^his correction term is both 
accurate and quickly computable. 

I'he vertical component of gravity from a curved triangular 
lamina segment is given in spherical coordinates by (see Figure 
5a) 

Ag - AZkp jj “1^ dA 


ff Af’sin fl(Z - cos tf) .. 

“ AZkp 1 1 ^ — i d9 d^ 


where 


k universal gravity constant: 
p density: 

SZ thickness of polygon: 

R radius of body: 

t distance from field point to edge of laminae: 

/, height of field pviinl above surface, equal loZ - R. 
a angle to vertical: 

0, d spherical coordinates. 


As R 


Ag — ♦ AZkp 




d^ 


( 2 ) 


the flat expression. 

Our task. then, is to compute the integral in (I) along 
each edge of the polygon. Unfortunately, unlike the ex- 
pression for a flat body [Talwani and Ewing, I960J. the 
expression for a curved body is not analytical! grable 
along an edge, nor is it easy to approximate by using scries 
or other integration techniques. However, the integral of 
the dilTcrcncc between the flat und the curved integrands can 
be easily approximated along an edge, und this diO'crence 
added to the exact flat expression results in a good approx- 
imation to the curved integral. 

To do this, we will divide each lamina segment into four 
equally spaced regions (in ^) (sec Figure ib) and apply a 
parabolic approximation to the integral using the resulting five 
points. I ct /(« be the distance from the nth point on the lamina 
to the vertical projection of the field point on the surface (see 
Figure 5t ). Let the dilVercnce between the itat and the curved 
expressions be 


D, I — 





(3) 


where r,' (/»’ ■* 4,’)* ’ is the flat distance from the held 

point to the nth edge point and r, j/f' < /' - 2RZ. ci»s 
{AkIR)V ' IS the corresponding curved disi.incc. Ilic integral 
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Fif. 5. Geometrical dcmenlt involved in calculaiinf the vertical 
pravnation attraction ol a ipherical triaiifular cap. 

of the dillcrence is then given b>' u parabolic approsimaiion 
for live equally spaced points; 

i - (0 96X90, + 280, + 220, + 2#0, 4 90,) (4) 

The 10011131) from each lamina segment is then 

where this integral is computed eaactly. 

It has been found, by using direct numerical integration, 
that fur must contours the accuracy of using an approximation 
at only five points is better than 0.3%. However, if the field 
point is above and near a contour {A^ -• 0) or if an edge is 
extremely skewed (0 • 0), these errors become significant. 
Fortunately, in the case in which Am is small enough to alTect 
these errors, the dilTerence between the curved and the flat 
values is less than 1% v '.he approximation errors are of the 

Older ol 50% and is less .■ 0.1% when the errors are as high 

as 100%. Thus the net result is an effective error of less than 
0.5%. Also, us 0 — 0, the effective area and hence the gravita- 
tional contribution approach 0. and so this source of error also 
becomes negligible. 

The testing wav done for a contour JOO k*n long. In practice, 
contours which the lield points will be near arc generally much 
shorter than this value, the errors thus being reduced to ihe 
accuracy of the machine. One could use seven or 10 points 
instead of live in the parabolic approximation, at small add!- 
tioiis ol lime, but exhaustive testing has shown this to be 
unwariaiilcd. 

Afknimlfjgnii'nit \ l.irgc numhei of dedicated people are revpon- 
mMc lur the ahitii) to uhlain information on the fravitv lield of the 
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Negative gravity anomalies on the moon 

M p . . - 

. , .wliv Carl Bowin 

Wood* Hole Occanosraphk IntliluUon,* Wood* Kol«, MBMachuMlli0234) 

Ab*lrKt>-Nef«liv« tnvily anoiratics on th« moon cm be catcjioriccd by whether they Jiow si 
correspondence to lunsu' topoyrephy or appear to be unrelated to surface topography. Clear examples 
of the former we the neyative anomalies from Doppler residual data that occur over craters such as 
Copernicus, Hipparchus, Ptolcmacus, Theophilus, and Lanjoeniis. These anomalies appear to be due 
to mass dcftciencies caused by the craterinf process, in larpe part probably by ejection of m:itcrial 
(tom the crater. Some anomalies on the fw tide deicrmir.cd from an analysis of orbital rate* (Ferrari. 
1975) may have a correspondence with topography, but others do not. Irregularities in the thickness of 
the lunar crust, not related to compensation for present surface topography, are inferred to be a likely 
source (or the broad Iwpe gruvity anomalies unrelated to surface topography. Localized large negative 
anomalies adjacent to mascont, such at these between Mare Sercnilatis and Mare Imbrium or Mare 
Critium, may be the result of localized depression of lunar crust adjacent to mate b.iiins during the 
cwly stages of mare filling. Structures on the moon having a half'Wavelcngth of COO km or lest, and 
large negative or positive gravity anomalies we not in isostatic equilibrium. M.iny of these features 
have mass loadings of about 1000 kg/cm* which can be statically sustained on the moon, even though 
such isostatic ma«s anomalies cannot be so maintained on the earth. 


Introduction 

Although thr large positivb lunar gravity anomalieg (mascons) have at- 
tracted the most attention, negative gravity anomalies (free-air anoinaties) have 
been recognized for the lunar surface since the first analysis of the gravity field of 
the moon (Muller and Sjogren, 1968) based on Lunar Orbiter data. The first 
free-air gravity anomaly map of the lunar front side, presented by Muller and 
Sjogren (1968), revealed minimum negative r.rR^ity anomalies of ’-lOOmgul 
northwest of Ptolemaeus. two lows of ~ 80 ingal between Mare Screntitatis and 
Mare Nectaris, a low of -70 mgal north of Serenitatis, two lows of - 60 mgal due 
west of Ptolemaeus, a low of -60 mgal northwest of Mare Imbrium. and several 
lows of lesser magnitude. The maKons shown on the same map are -f 60 mgal for 
Mare Humorum, + 100 mgal for Nectaris and Crisium, 4- 180 mgal for Serenitatis. 
and 4- 230 mgal for Imbrium. Two highs of 4 70 mgal apparently unrelated to 
circular maria were also identified, one southeast of Eratosthenes and the other 
farther southeast, but northwest of Sinus Medii. *l'he closest approach of the 
Lunar Orbiter missions to the lunar surface was 100 km, and the anomalies were 
normalized to an altitude of 100 km under the assumption that mass variations 
were located at a depth of SO km below the surface. Information from thc^Afollo 
missions has confirmed the existence of mascons, and led to revision of some of 
the other features. The only clear relationship of the orbiter gravity anomalies to 
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Fig. 1. Frce-air gravity anomalies from Apollos 15 and 17. Frofiles of total acceleration 
(estimate of free-air anomaly) obtained by processing Doppler data using extendetl 
sequential estimation techniques by W. Wollenhaupt at Johnson Spacecraft Center. 
These profiles are presented to show the attenuation of the r‘.f>ative anomalies with 
increasing height. Profile at 300-km height from Apollo 15 Kev I, ;it lOO-km height from 
Apollo 15 Rev. 19, and at 20-km height from Apollo 17 Rev. 3. 


topography is the anomalous association of the mascon anomalies with depres- 
sions in the surface at the circular maria. It was not until the later Apollo missions 
that gravity anomalies were determined from orbits which approached to within 
about IS km of the surface. In some of those revolutions of the Apollo cpmmand 
ship module (CSM), negative gravity anomalies were observed to correlate with 
the locations of craters having diameters of less than 100 km. Thus, both gravity 
anomalies that correlate with topography, cither directly or inversely, and those 
that do not are present on the moon. ^ 

Iligh-resolution, global coverage of the lunar gravity field is not yet available. 
Interpretation rests upon the extremely limited surface measurements, the 
incomplete coverage of the moon’s near side by high- and moderate-resolution 
gravity data obtained from analysis of Doppler residucals for individual spacecraft 
orbits (e.g. Fig. t), and the lower resolution, more complete global coverage (Fig. 
2) obtained from analysis of orbital rales (I'crrari, 1974, 1975). A similar global 
• analysis also has been prepared by Ananda (1975). 
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Fig. 2. Regional frec-air gravity anomalies of the moon and eArth. A is a fifteenth qrxler 
and degree spherical harmonic free-air anomaly map at 100-km height of the liinar 
surface from Ferrari (1973). B is a regional free-air anomaly map of the earth at the 
surface from Smithsonian Astrophysical Observatory (SAO) Standard Earth II spherical 
harmonic coefficients referenced to the International Gravity Formula 1930: approxi- 
mately a sixteenth order and degree map. Contours in milligals. 
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Negative Anomalies Related to Topography 

Negative anomalies over unfilled craters have been observed in the gravity 
data obtained from tracking the low-altitude perinpsis (less than SO km) orbits of 
Apollo missions. These generally occurred during the early revolutions (revs.) of 
the Apollo missions (e.g. Fig. 1; 20-km height), prior to separation of the LEM, 
and more continuously in the case of the subsatellites of the Apollo 15 and 16 
missions. Examples are the — S7-mga1 anomaly determined at about 20-km altitude 
above the edge of Copernicus Crater (Sjogren et al., 1973, 1974), — 67mgal at 
about 16-km altitude nearly above the center of Hipparchus (Sjorgen et al, 1972), 

— 9S-mgal anomaly determined at about 30-krn altitude nearly over the center of 
Ptolemaeus (Sjogren, -1974), — 120mgal at about 16 km over the center of 
Theophilus (Sjogren et al, 1972), and the - 127-mgal anomaly determined at about 
2S km above the edge of Langrenus (Sjogren, 1974). 

These anomalies appear to be clearly due to mass deficiencies caused by the 
cratering process, in large part probably by ejection of material from the crater. 
Disk-shaped models representing the deficiency of crater depressions have been 
successful in reproducing the general magnitude of the observed gravity 
anomalies (Sjogren et al, 1972; Scott, 1974). Sjogren et al (1974) found that their 
best fit to the magnitude and gradient of the Copernicus anomaly was obtained 
when disk masses were included to simulate the mass of the ciater rim beneath the 
orbit trajectory. 

Negative anomalies of much broader extent than typical individual craters are 
common on the moon, both on the near side and on the far side (Fig. 2A). Those 
on the near side show little or no correspondence with topographic depressions. 
Some of those on the far side have been inferred by Ferrari (1974, 1975) to be 
associated with ringed biisins such as Mendeleev, Moscoviense, Korolov, and 
Apollo, and he associates other ringed basins such as HertZM ’‘ung. Ingenii, 
Gagarin, and Tsiolkovsky with relative lows in regions of gravity highs. The 
correlation of Korolov with a -110-mgal minimum and that of Apollo with a 

— 166-mgal minimum are the best. However, the size of the ringed basins as 
presently known; the lack of similar anomalies over other ringed basins such as 
Leibnitz and Milne; the apparent lack of ringed-basin structures correlative with 
the - 112-mgal (15N,255E), -93-mgal (15N,200E), - 130-mgal (45S,165E), and 

— 145-mgal (45N, I95E) minima; and the different local gravity relief for similar 
sized features such as Mendeleev, Gagarin and Ingenii, and Hertzsprung and ' 
Apollo suggest that the broad negative anomalies of the back side, also, are not ' 
simply a result of topographic depressions. More refined processing of the orbital 
rate data may contribute to more definitive knowledge of the location and 
magnitude of the regional gravity anomalies. 


NiiOATIVE ANOMAt,IES NOT RELATED TO TOPOGRAPHY 

The cause of these negative gravity anomalies is not obvious. They do not 
appear to be related to topographic depressions, and commonly occur over 
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highland topography. The most prominent negative anomalies encountered in the 
Doppler residual acceleration data are those between the Mare Imbrium and Mare 
Serenitatis mascon anomalies, and between the Mare Serenitatis and Mare 
Crisium mascons (Fig. >l);.. Respectively, i. ese anomalies are about -20 and 
— 5 mgal at 300*km height, about — SO and — 60 mgal at 100-km height, and about 
-250 and —225 mgal at 20-km hei^.t '^»o«r«*n et ai, 1974; Bowin ct al., in press; 
Muller et al., 1974). Figure ! shows a sampling of the Doppler residual data for 
those heights but do not necessarily show the most pronounced values attained. 
Bowin et al. (in press) demonstrated that the gravity low to the east of Mare 
Serenitatis does not continue along the southern margin of the basin. Thus these 
negative anomalies apparently do not form a ring :uound Mare Serenitatis and 
very probably neither do they ring the other circulsr mare basins. Furthermore, 
they are too negative to be the result of an edge effect from a deep mass deficiency 
that might result if a root of crustal material existed beneath the mare basins to 
provide isostatic compensation, as suggested by Hulme (1972). Instead of a crustal 
root, a mantle dome appears to be more likely beneath the mascons. The more 
recent proposals for the description of mascons have been thick mare fillings 
(Phillips et al., 1972), mantle plugs (Wise and Yates, 1970), and a combination of 
mantle dome and thin mare filling (Bowin et al., in press). The combination 
explanation appears to best fit the existing gravity and geologic information. The 
mascons require a mass distribution that cannot be in isostatic equilibrium, and 
hence they indicate a strong lunar crust. 

Possible explanations for the negative (and positive) gravity anomalies that are 
not associated with surface topography could be (1) lateral variations in crustal 
density, (2) lateral variations in the density of the lunar mantle, or (3) irregularities 
in the thickness of the lunar crust not related to compensation for surface 
topography. Although conceivable, lateral density variations of the magnitude 
required to account for the broad, large negative gravity anomalies (or positive 
anomalies not associated with circular maria) are not considered as likely as 
irregularities in the thickness of the lunar crust. 

The cause of the suggested uncompensated irregularities in the thickness of 
lunar crust is speculative. Roughly, about 1-km variation of surface height, if 
uncompensated, would produce a free-air gravity anomaly variation of about 
120 mgal. Assuming a lunar crust-mantle density contrast of 0.5 g/cm*, then about 
5.8 km of subcrustal variation in depth to mantle would produce about the same 
gravity anomaly. Erosion of a kilometer or two of surface material, as the result of 
impacts when the crust was rigid enough to resist isostatic adjustments of the 
affected areas, could possibly produce some of the observed negativej gravity 
anomalies. If we imagine that such erosion and attendant redeposition might tend 
to smooth overall the lunar topography, then the gravitational contribution of the 
subcrustal depth variations might become evident in the external field, and 
produce some of the broad anomalies observed. Both negative and positive 
regional gravity anomalies that would not correlate with surface topography could 
be produced in this manner. 

To explain the large negative anomalies adjacent to mare basins. (I'ig. I), the 
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following scenario i‘ offered. An impact created a mare basin, beneath which a 
mantle dome rose to provide isostatic compensation for the deficiency of the 
basin above. It is also inferred that the highland rim surrounding the basin may 
have been higher than it is today, and that it was isostatically compensated by a 
crustal root into the mantle beneath. The root was probably deeper in some places 
than in others around the periphery of the mare basin, and these places may be the 
sites from which more mantle material flowed during the rise of the mantle dome. 
It is possible that these places have greater crustal thickness because more ejected 
mcaterial was deposited there. Evidence from crater counts and radioactive age 
dating of returned samples suggests that considerable time (perhaps about 
200xl0*yr) elapsed between the formation of the. large mare basins and their 
filling with mare basalts (Baldwin, 1963, 1971; Huneke et al., 1973). During this 
time interval, the rigidity of the lunar crust is inferred to have increased (Phillips 
ef at., 1972; Baldwin, 1971; Miiehlberger, 1974). The initial outpourings of mare 
fill, which has a density about O.S g/cm* greater than that of lunar crust (Bowin et 
at., in press; Talwani et at., 1973), may have been compensated by subsidence of 
the mare basin.' Since to a first approximation the mantle dome has the same 
density as the mare fill, it would not rise further to provide isostatic compensation 
for the subsidence of the basin floor due to the loading of the fill. The top of the 
dome would sink the same amount as did the basin floor. Faults at the edge of the 
mare fill could have isolated the adjacent crust from the subsidence taking place 
within the basin. However, as the rigidity of the lunar crust increased, local 
isostatic compensation could no longer be complete, and further subsidence of the 
basin, due to additional filling, would begin to depress portions of the crust 
surrounding the basin. Continued faulting at other portions of the basin margin 
may have kept their adjacent rims from being depressed. The inner edges of the 
shelf areas of mare basins may indicate the location of some of the sites of such 
concentric faulting that occurred during mare filling. Perhaps the thickness of the 
crust beneath the rim partially determined the degree to which different portions 
of the rim were depressed. Depression of crust that formerly was in isostatic 
equilibrium, of course, produces a mass deficiency and hence negative gravity 
anomalies. The more localized the area of depression, the greater will be the 
gravity gradients, and the greater the depression, the greater the negative 
anomalies. Eventually, the lunar crust became rigid enough that it could support 
the last fill without subsidence of the basin occurring, and, thus the excess mass of 
the fill together with the mass of the mantle dome produces the mascon gravity 
anomalies. • \ 


ISOSTASY 

Kaula (1969) argued that the moon is closer to hydrostatic equilibrium than the 
earth because the lower degrees of the power spectrum of the lunar gravity field 
have lower absolute magnitude than that predicted from the earth’s gravity field 
under the equal-stress assumption, and because the excess mass of the largest 
lunar mascon is smaller than that of a score of large positive anomalies on the 
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earth. Although the power in the low degrees is lower for the moon than for the 
earth, the power in the higher degrees (7-13) was determined to be higher for the 
moon. This relation for the higher degrees was taken by Kaula (1969) to be a 
consequence of errors and distortion in the analysis of the data. Figure 2, 
however, suggests that that relation is largely real and hence has physical 
significance. It is of course now desirable that the power spectrum be redeter- 
mined utilizing more recent data. O’Keefe (1968) calculated Bouguer anomalies 
for the lunar near side and concluded that, to a rough first approximation, the 
moon is in isostatic equilibrium. Baldwin (1971) seems to argue that the moon is 
nearly in isostatic equilibrium, including the mascons, because it approaches the 
condition of isostatic equilibrium as closely as the threshold strength of the lunar 
materials permits. • 

Figure 2A and the power spectrum data presented by Kaula (1969; Table 3) 
suggest, I claim, that structures on the moon having half-wavelengths of approxi- 
mately SOO km or less are commonly not in isostatic equilibrium. These features of 
the gravity field represent mass loadings or deficiencies approaching about 
1000 kg/cm*, which appear to be statically maintained by the strength of the outer 
part of the moon. Obviously, the lower the gravitational attraction a planetary 
body has,' the greater a mass per unit area loading it can sustain for material with 
equivalent strength. Thus, isostatic mass anomalies of 1000 kg/cm* can be 
statically sustained by earth-like materials on the moon, or Phobos, or Icarus, for 
example, but not on a body with the gravitational attraction of the sun where 
statically maintained isostatic anomalies of only a few kilograms per square 
centimeter would be expected. Neither, apparently, can so large a load be 
sustained on the earth. The largest isostatic mass anomalies on the earth appear to 
be about 1000 kg/cm* (Hess, 1935) and are associated with island arcs. The gravity 
anomalies associated with these featur<>s are not statically maintained but 
apparently decay rapidly upon cessation of the dynamic conditions that produce 
the island arcs and adjacent deep sea trenches. Furthermore, the half-wavelength 
of arc-trench structures is only about 200 km at the maximum. 

The lunar gravity anomalies shown in Fig. 2A are much greater in magnitude 
than earth anomalies averaged over equivalent spacial dimensions. The fifteenth 
order and degree lunar anomaly map (Fig. 2A) has local gravity relief commonly 
of over 300 mgal. Fifteenth order on the moon is approximately equivalent to a 
three degree square average on the earth. In the Caribbean region,* which has 
perhaps the greatest local gravity anomaly relief on the earth, only a little over 
100-mgal relief is attained in a three degree square average (Bowin, in press). 

I conclude that many structures on the moon having a half-wavelength of 
800 km or less are not in isostatic equilibrium, and that mass loadings of about 
1000 kg/cm* cannot be sustained statically on the earth as they are on the^moon. 
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■ 'CbMPARiSON OF GRAVITY ANOMALIES FOR EARTH, MARS, AND MOON 
Carl Bowln, Woods Hole Oceanographic Institution 
Woods Hole, Massachusetts 02543 

It is only within the last eight years that man has been able to com- 
pare gravity anomalies with topography on planetary bodies other than his 
own. Topography shows the distribution of surface irregularities, and gravity 
anomalies show the distribution of mass anomalies. The coherence, or lack 
of it, between these two parameters provides significant information on 
crustal structure and its evolution. 

On the earth at small scales (less than 100 to 300 km half wavelength) 
there is generally a positive correlation between free-air gravity anomalies 
and topography. However, towards the larger dimensioned features, the free- 
air anomalies are commonly reduced in amplitude indicating partial to com- 
plete compensation for the mass of the topography. From gravity anomalies 
and topography the thickness of the elastic part of the lithosphere may be 
estimated (1). These studies suggest a thickness for oceanic lithosphere of 
about 10 km near spreading ridges, to values near 30 km for 60 million year 
old lithosphere to about 41 km in older lithosphere at deep sea trenches. 
Beneath continents the elastic part of the lithosphere may be about 60 km 
thick (2). Mass loadings of up to about 400 kg/cm^ can be statically main- 
tained on the earth as with the example of the Ninety East Ridge in the 
Indian Ocean (3). The greatest mass anomalies are associated with trenches 
where mass deficiencies approaching 1,000 kg/cm^ occur (4). As pointed out 
previously (5), the gravity anomalies associated with trenches are not 
statically maintained, but decay rapidly upon cessation of the dynamic con- 
ditions that produce island arcs and adjacent deep sea trenches. On the 
earth, the largest gravity anomalies are associated with features that have 
half wavelengths of less than 200 km and that are dynamically maintained by 
the plate tectonic processes active today. 

At large scales (greater than 100 to 300 km half wavelength) there is 
generally a lack of coherence between free-air gravity anomalies and topo- 
graphy, and the anomalies over the vast majority of the surface are small in 
magnitude (less than 30 mgal for spherical harmonic models of degree and 
order 16)(6), indicating rather close adjustment to isostatic equilibrium. 

Mars has a radius of only about half that of the earth, and a surface 
gravitational attraction only 387. as great. However, Mars has very large 
regional gravity anomalies. Mariner 9 Doppler data provides information on 
the gravity field variations to the resolution of an eight or nine degree and 
order model (7, 8, 9, 10). Measured free-air anomalies range from -300 mgal 
to +500 mgal (11). The large martian gravity anomalies correlate well with 
regional topographic variations, the high positive values occurring over the 
Tharsis Plateau region, and the negative values over the Chryse and amazonls 
lowlands. The anomalies, however, ore only about four-tenths ns great as 
they would be if the topography were not partially compensated (11). The 
remainder of the martian topography, at least at the resolution of 8th order, 
is considerably nearer to being in Isostatlc compensation, although anomalies 
50 to 100 mgal from zero are common. The isostatic deviation map prepared by 
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Phillips and Saunders (11) Ividlcates a 1,200 kg/cm^ mass anomaly associated 
[w^th the Tharsls Plateau region. The mass anomalies associated with the 
gla..* shield volcanoes are not yet resolvable. 

The moon's gravity anomalies Indicate that structures having: half-wave- 
lengths of approximately 800 km or less are commonly not in isostatlc equlll- 
jbrlum (5). Globally Its gravity anomalies have been determined to a spherical 
harmonic model representation of fifteenth order and degree (12). Free-air 
gravity anomalies 't 100 km altitude range from -166 to +260 mgal. Both these 
extremes are locatea the lunar farslde where control Is the weakest. The 

range of frontside anou.. Hes is from -137 to +239 mgal. At 100 km altitude, a 
fifteenth order model for the earth shows a range of free-air anomalies from 
-0.1 to +0.1 mgal. On the lunar nearside the large positive anomalies are all 
associated with mascons and indicate mass anomalies of about 800 kg/.cm^( 13,14), 
The available data for three planetary bodies indicates the following 
'comparison. The magnitude of the maximum gravity anomalies and the mass per 
unit area anomalies differ by less than a factor of two amongst the three 
bodies. The areal size of the maximum anomalies, however, differ much more 
^rkedly. On the earth they are restricted to half wavelengths of less than 
200 km. On the moon they extend up to half wavelengths of 800 km. But on Mars 
they extend to over a 4,000 km half wavelength.' On the earth, the largest 
gravity anomalies require dynamic processes to maintain them. Maximum static 
jloads that can be sustained appear to be on the order of 500 kg/cm^, and this 
jlow value is presumably due to the thinness of the rigid lithosphere. For the 
'earth and Mars, the large gravity anomalies generally are correlated with 
topography. The mascons on' the moon, however, are inversely correlated with 
topography, the large positive anomalies occurring over mare basins on the 
^earside. The large anomalies determined for the farside have been inferred 
jto be positively correlated with regional topography irregularities (12) but 
this generality has been questioned (5). Visual inspection of the topographic 
land gravity Lunar Orbital Data Maps recently published (13, 14), allowing for 
jthe different resolutions, does not support a strong correlation of topography 
and gravity on the farside. The planned Lunar Polar Orbiter mission should 
resolve the matter. Each of the three bodies appears to have experienced a 
^>eriod during which its crust was largely in isostatic equilibrium. The moon 
left this stage about 3 b.y. (17), Mars developed the Tharsis Plateau perhaps 
within the last 1 b.y. (11, 18), and the earth's crust remains preponderantly 
lin isostatic equilibrium. Although the mascon mass anomalies lie within 100 
km of the lunar surface (13,14), it is no doubt a thick rigid lithosphere that 
Isupports them. Although the Tharsis mass anomaly on Mars is greater than the 
ilunar mascons, this does not necessarily require a thicker rigid lithosphere 
on Mars. It seems more likely that on the moon there was simply not suffi- 
cient extrusion of dense volcanic material to create larger mascons. 

The gravity data is compatible with the moon. Mars, and the earth being 
examples of an evolutionary sequence. The moon became frozen at the youngest 
age, perhaps while bombardment was still at n moderately high flux.. Mars de- 
veloped greater volcanic activity, supplying greater volume to the outer part 
of the crust than did the moon, following development of a thick lithosphere. 
The Tharsi.s Plateau region may be a gigantic lacolith. The shield volcanoes 
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of Mars may be terminal volcanic outbursts. The earth's heat engine appears to 
have had the greatest sustaining power and plate tectonic activity developed, 
and a thin lithosphere persists. 

Perhaps Mars offers a clue to the origin of the asteroid belt. I suggest 
that a planetary body that formerly existed there, in evolving through the 
volcanic stage developed a surflclal mass anomaly significantly larger than, 
that observed in the Tharsls Plateau. Sjogren et al. (7) found that a single 
anomalous mass about 0.8 x 10"^ ^ars beneath the Tharsls region at a 

radial distance of 0.571 %iars accounted for the main features of the ano- 
malous gravity field. We may take this as a reasonable estimate of the mass 
of the Tharsls Plateau anomaly, but its actual depth would be much shallower 
because of Its broad special dimensions. A larger mass anomaly in the outer 
portion of a planetary body might lead to instabilities of its motion causing 
its break-up. 
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